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Abstract 

We review the formalism by which the tunnelling probability of an unstable ground state can be computed in quantum field theory, 
with special reference to the Standard Model of electroweak interactions. We describe in some detail the approximations implicitly 
adopted in such calculation. Particular attention is devoted to the role of scale invariance, and to the different implications of scale- 
invariance violations due to quantum effects and possible new degrees of freedom. We show that new interactions characterized 
by a new energy scale, close to the Planck mass, do not invalidate the main conclusions about the stability of the Standard Model 
ground state derived in absence of such terms. 
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1. Introduction 

In recent years there has been considerable interest in the 
problem of the stability of the Standard Model (SM) ground 
state. Due to the sizable negative contribution to the /? func¬ 
tion of the Higgs self-coupling induced by top-quark loops, 
the usual electroweak vacuum |0), characterized by (0|/!|0) = 
V =» 246 GeV, may not be the absolute minimum of the scalar 
potential. In this case the true minimum of the theory is lo¬ 
cated at much larger energy scales and, in the absence of New 
Physics (NP) modifying the effective Higgs potential, the elec¬ 
troweak vacuum is unstable. The parameter space of the SM 
(with particular reference to the top-quark mass, m,, and the 
Higgs-boson mass, m*, which are the most relevant parame¬ 
ters) is thus naturally divided into three regions: stability, insta¬ 
bility, and metastability. The stability region is the one where 
the electroweak vacuum is the absolute minimum of the poten¬ 
tial. The instability and metastability regions are those where a 
new deeper minimum exists, with the metastability region being 
characterized by a lifetime of the unstable electroweak vacuum 
larger than the age of the Universe. More precisely, the insta- 
bility/mestastability boundary is determined by the decay prob¬ 
ability of the electroweak vacuum under quantum tunnelling, 
that sets a model-independent upper bound on the lifetime of 
the unstable vacuum irrespective of the thermal history of the 
Universe. 

A precise determination of the boundaries of these three re¬ 
gions has recently been presented in Refs. |[T1 |3 (see also 
Eiiiiia). As a result of these recent analyses, the present exper¬ 
imental values of m/, and m, lie in the metastability region of the 
SM parameter space. This finding holds only within the SM; 
however, it has an important consequence for beyond-the-SM 
searches: it implies that there is no need to invoke the presence 
of NP in order to stabilize the SM electroweak vacuum. 

The validity of the analysis in Refs. |[T1 |2l and the corre¬ 
sponding conclusions has been questioned in a series of recent 


papers EllTlElEl. There, it has been shown that non-standard 
physics modifying the shape of the Higgs potential at energy 
scales of the order of the Planck mass can sizeably affect the 
tunnelling rate of the electroweak vacuum. This observation is 
correct. However, as we discuss in the following, it does not 
invalidate the interest and the main conclusions of the analyses 
based on the SM potential. 

In this paper we present a critical re-analysis of the problem 
of the SM vacuum stability. Our purpose is to clarify the as¬ 
sumptions and approximations employed in the evaluation of 
the SM tunnelling rate, with particular attention to those that 
have been often overlooked or implicitly adopted in the exist¬ 
ing literature. 


2. The tunnelling rate within the SM 


We begin by reviewing the standard formalism, originally 
worked out by Coleman and Callan IfTOlfTTI . which allows one 
to compute the probability per unit time (or, equivalently, the 
lifetime) of a false ground state to a true ground state in quan¬ 
tum field theory. In the semiclassical approximation, the de¬ 
cay probability per unit time of the electroweak ground state is 
given by lfT2l 

r (1) 

where tu is the age of the Universe, 


5[/r] = 



df,hd,,h + V(h) 


( 2 ) 


is the euclidean action of the theory, computed for a specific 
solution h of the euclidean field equation for the scalar field 
which is usually called the bounce, and R is a dimensional fac¬ 
tor associated with the size of the bounce. The bounce field 
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configuration is such that it is equal to the false vacuum conhg- 
uration h = v at inhnite euclidean time t, and completes barrier 
penetration at t = 0. 

A known result, conjectured by Coleman and subsequently 
proved by Coleman himself, Glaser and Martin ca, guarantees 
that the bounce solution of minimum action is invariant under 
four-dimensional rotations in euclidean spacetime, that is 

h - h(ry, (3) 


Hence, 


lim h(r) - v. 

r-^oc 


(4) 


By requiring that the solution is non-singular at the origin, we 
also have 


dh{r) 


dr 


= 0 . 


r=0 


(5) 


Eq. (jT]) gives the leading contribution to the tunnelling rate in 
the semiclassical limit, that is, it only includes exponentially 
enhanced terms in the limit —» 0. In particular, the over¬ 
all normalization can only be determined by including the first 
quantum corrections im. 

We observe that the bounce is the unique solution of a suit¬ 
ably defined Cauchy problem. Indeed, under the assumption 
that the bounce is (9(4)-invariant, the field equation for the 
bounce is 

h"{r)^-\'{r) + V'(h), (6) 

r 


where V(/i) is the scalar potential, and primes denote differen¬ 
tiation with respect to the functional argument. A set of initial 
conditions 

d(ro) = ko; //(ro)^hi, (7) 


given at any hnite value ro + 0, dehnes a Cauchy problem, 
which has a unique solution in a neighbourhood of ro, as a con¬ 
sequence of known results in real analysis. Since r = 0 is the 
only singular point in the r.h.s. of Eq. ro can be taken to 
be arbitrarily large. However, the unique solution is not nec¬ 
essarily well dehned at r = 0, because of the singularity in the 
r.h.s. ofEq. nor at r ^ oo, because the existence and unicity 
theorem has a local meaning. Thus, for a generic choice of the 
initial conditions at r = ro, the boundary conditions Eqs. ®, 
are not necessarily fulhlled. Conversely, the requirement that 
Eqs. are fulhlled by the bounce selects a set of allowed 
initial conditions at an intermediate point ro. Depending on the 
shape of the potential, it is possible to have multiple bounce so¬ 
lutions (with different initial values !iq and hi) which satisfy the 
boundary conditions in Eqs. i. 


2.7. Decay of the SM vacuum in the semiclassical approxima¬ 
tion 

The case of the pure SM is especially interesting. In this 
case, the value of the true vacuum is typically very large with 
respect to the electroweak scale v a; 246 GeV, so one usually 
takes V = 0. Thus 

lim h(r) = 0. (8) 


The validity of this approximation is discussed below in 


Sect. 2.5 The scalar potential in the unstable region is there¬ 
fore 


V(h) = ^Ah\ 


(9) 


where we neglect the logarithmic running of A and take it as a 


negative constant. This approximation is reviewed in Sect. 2.2 
Then Eq. (j^ takes the form 


h''{r) H— h'(r) - Ah^(r). 
r 


( 10 ) 


Eq. (10 1 is invariant under scale transformations; if h(r) is a 


solution, then 

ha(r) — ah{ar) (11) 

is also a solution, for any choice of the scale factor a. Indeed 


Kir) + -h'^{r) = a^ 


h"{ar) H- h'(ar) 

ar 

A 1.3 


= Aa^hfar) = T/r^(r). (12) 


Obviously, the scaled solution ha{r) has the same limiting be¬ 
haviors Eqs. w as the original one, but different initial con¬ 
ditions at r = r(j. Otherwise stated, the boundary conditions do 
not fix the overall normalization of the bounce solution. 

As is well known, a solution of Eq. ( [T0) l with the boundary 
conditions 03 

na 


given by the Eubini-Lipatov instanton lfT4l 


hir) = 


8 R 


(13) 


|4| 7^2 -H r2 ’ 

for any value of R and 4 < 0 (note that h(R) — h(0)l2: this 
will be our dehnition of the size of the bounce throughout the 
paper). It should be clear from the above discussion that the 
presence of the arbitrary parameter R is just a reflection of the 
scale invariance of the equation, and not a signal of non-unicity 


of the solution; indeed, the scaling defined in Eq. (Ill amounts 


to replacing R with R/a in Eq. GD- h(r) is the unique solution 
of Eq. 0 with initial conditions in ro 


ho — 


8 R 
7^2 


h] = 


8 2Rro 
| 4 | ( 7^2 -H ^ 2)2 ’ 


(14) 


which obey the constraints (^i. 

The SM bounce Eq. ( [T3] l can be found by the following pro¬ 
cedure (see e.g. Ref. IfThI for an alternative derivation). Let us 
assume that a solution of Eq. ( fTO] ! exists, with a Taylor expan¬ 
sion around r = 0; 


h(r) = ^ Akd 

k=0 

with Ao > 0. Eq. ( [T0| i takes the form 
3A 

- -I- 

k=0 

It follows that 


(15) 


k k - 


- + + 2)(k + 4)A,+2r'^ = 4 AiAjAk-. 

^=0 1=0 j=0 


A, = 0. 


i-j- 

(16) 

(17) 
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Thus, the condition that the first derivative of the bounce at r = 
0 vanishes is a consequence of the assumption Eq. ( [T5| ). The 
remaining coefficients are given by the recurrence relation 

" {k + 2){k + A) § 2 (18) 

The coefficients Ak with k odd are zero; indeed, when k is odd, 
one out of the three summation indices i, j, k - i - j is also odd 
(possibly all of them). Hence Ai = 0 implies A 3 = 0, and so on. 
Thus, we may rewrite Eq. ( [T5| ) as 

00 

h(r) = ^ ak = A 2 k, (19) 

A:=0 


and the recurrence relation for the coefficients becomes 

= 8 WTWTT) § § ( 20 ) 

The coefficients at are determined by the single number ao, the 
value of the bounce at the origin. We now show that 


The value of the euclidean action S [h] of the bounce solu¬ 
tions Eq. (fTSll is 

87r2 

sm---, (26) 


independently of the value of R. This is not surprising, because 
the action is dimensionless in natural units, and no dimension- 
full scale parameter is available other than R. Hence, there is no 
way to single out one preferred value of R at the semiclassical 
level. However, R is related to the value of the bounce at r = 0: 


and since the bounce solution only exists for d < 0 , we expect 
that 

1 > A/, (28) 

K 

where A/ 10'° GeV is the energy scale at which the running 
coupling /l(ju) becomes negative. This is an a posteriori confir¬ 
mation that neglecting the electroweak scale, of order 10^ GeV, 
with respect to the size of the bounce is indeed a reliable ap¬ 
proximation. 




2j+l 


The proof is by induction. Eor k - 0 Eq. (|20|l gives 


( 21 ) 


Cl\ 


8 


( 22 ) 


We now assume that Eq. (21 1 holds for 0 < j < k. Then 
A 2 


aic+\ 


k k~i 

y, y, at 


8(k+l)(k + 2)^^ 

k k-i 


Q jQ.k~i— j 


8 (^+l)(^ + 2 ) 




,2k+3 




^2k+'i 


(23) 


which is what we set out to prove. The Taylor expansion in 
Eq. ([T9)| can now be summed. We find 


h(r) - ao 


go 


JlkJlk 


ao 


l-^alr^' 


The series has convergence radius 




|d| ao 


(24) 


(25) 


but if A is negative the sum can be analytically continued to the 
whole positive real axis, and vanishes as r — » oo. Finally, we 
note that for d < 0 the bounce in Eq. ( |24] i coincides with the 
solution given in Eq. (©, with R as in Eq. ( [25] l. The above 
construction shows that, given the value of the bounce at the 
origin and the requirement of regularity on the range 0 < r < oo, 
the solution is unique. 


2.2. Violation of scale invariance through radiative corrections 

The first quantum corrections, computed in Ref. HI, affect 
the semiclassical result in two respects: they fix the normal¬ 
ization in Eq. and they take into account the running of 
the Higgs coupling d[^ As a consequence, the tunnelling de¬ 
cay rate is dominated by the bounce with the maximum value 
of |d(ju)| ~ |d(l/RsM)l- For the central values of the SM pa¬ 
rameters the scale 1/Rsm turns out to be a couple of order of 
magnitudes below the Planck mass Mp = 1.22 x lO'® GeV. 

Schematically, one expands the euclidean action around the 
tree-level bounce solution 

S[h + h]^ S[h] + ^ J d^xS"[h]P, (29) 

and integrates over the fluctuations h: 

J =se-^W(det5"[/!])’^ (30) 

In the treatment of the functional determinant there are two 
main aspects which eventually lead to the appearance of the 
scale p ~ 1 /Rsm in the calculation; i) ultraviolet (UV) di¬ 
vergences of the non-zero modes of S [/i], responsible for the 
introduction of a renormalization scale p and ii) the treatment 
of the zero modes of S "[h], which have to be singled out and 
treated separately in order to avoid unphysical divergences. The 
existence of zero modes is simply a reflection of the fact that 
the classical action is invariant under a larger class of symme¬ 
tries which are broken by the explicit solution of the tree-level 
bounce (e.g. in the case of Eq. ([T3]l these are translations in 


* For some earlier works discussing the breaking of scale invariance due to 
the running of A in the context of the SM-vacuum tunnelling calculation see 
also Refs. iniEa. 
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0(4) and scale transformations). The symmetries of the action 
are hence restored only if one considers the family of all bounce 
solutions as a whole. This defines a measure of integration (the 
instanton measure m) over the the so-called collective 
coordinates, parametrizing the families of bounces with differ¬ 
ent sizes (R) and located anywhere in euclidean spacetime (xq) 


dfii, 


dRcfxQ 
R^ 


■ exp 




3|T(1//?)| 


-^S 


(31) 


where /l(l/f?) is the renormalized coupling at the scale IjR, 
AS - 0(^gf/A^ (with i running over the SM couplings) comes 
from the inclusion of the non-zero modes and contains finite 
terms plus logs which are minimized by the choice fi ~ IjR 
IfTSlI . The tunnelling probability is therefore obtained by inte¬ 
grating over the collective coordinates 


r , A cdR 


1 £/A(inst~Ty 1 ^exp 

AS 

3|4(1/R)| 


where the integral is extended to the range where /l( 1 /R) < 0, 
where a bounce exists. 

By expanding A( 1 /R) < 0 around its minimal value located 
at = Rsm 


A(l/R) ^ A(1 /Rsm) + ^T"(1//?sm)( 1//? - l//^SM)^ (33) 


we can approximate the integral by the method of the steepest 
descent 


4 3d2(i//;sM) 

^SM y 


exp 




3\A(1/Rsm)\ 


-AS 


^SM 


exp 




3\A(1/Rsm)\ 


(34) 


where in the last step we have neglected the subleading 
pre-exponential factors (including AS), and we have set 
/1"(1 //?sm) ~ dimensional grounds. 

The result Eq. ( [34| provides the standard approximation 
which is usually employed in order to estimate the lifetime of 
the electroweak vacuum. Nonetheless, it is interesting to com¬ 
pare it with a more direct calculation obtained via the renormal¬ 
ization group (RG) improved effective potential 


Reft = -^Aeff(h)h'^ 


^-A(h)h‘^ 


(35) 


where, as in Eq. we neglected the mass term and in the last 
step we have approximated the effective quartic coupling Teff 
with the MS renormalized coupling A evaluated at fi - h. The 
scale ambiguity of the tree-level bounce is now resolved by the 
leading order effective potential which takes into account the 
dominant radiative corrections. However, since it is not possi¬ 
ble to find an analytical solution for the bounce one has to resort 
to a numerical analysis 


plying our numerical 


;is (se^ 
set-uf[^ 


[Appendix A for details). By ap- 
we obtain the bounces displayed 


^In this work we refer to arXiv v4 of (2) for the central values of the SM 
couplings evaluated at the top-quark pole-mass scale. In particulai; we take 
mt= 173.34 GeV andm/, = 125.15GeV. 


in Eig.[^where, for illustrative purposes, we considered the case 
where the running of A is evaluated at one (It), two (2€) and 
three (3T) loops, respectively. The corresponding bounce ac- 



Figure 1: Bounce profiles for the potential in Eq. j35h where the running of A 
is evaluated at one (black), two (blue) and three (red) loops, respectively. The 
adimensional variables are normalized in terms of the Planck mass Mp. 


tions are found to be 

S[h,A(hy^] 

S[h,A(h)^^] 

S[h,A(h)^^] 


= 772.3, 

(36) 

= 1703.9, 

(37) 

= 1788.8. 

(38) 


On the other hand, we want to compare the latter case with the 
standard method (cf. Eq. p^ ) where A is evaluated at the scale 
1 /Rsm where [Sx = dtogfi ~ parameters, we 

l/R^^^ = 4.145 X 10^’ GeV; 4(1 /Rsm)'^ = -0.03409, 
1/R^^ = 2.768 X 10^’ GeV; 4(1 /Rsm)^^ = -0.01546, 
1/R^^ = 2.769 X 10^’ GeV; 4(1 /Rsm)^^ = -0.01473, 

which yields the bounces displayed in Eig. Correspondingly, 
the bounce actions defined by 


5 [/i, 4(1 /Rsm)] 


3|4(1/Rsm)I 


(39) 


have the following values 

5[/i,4(1/Rsm)''] =772.0, (40) 

5 [/i, 4(1 /Rsm)"^] = 1702.6, (41) 

5 [/i, 4(1 /Rsm)^3 = 1787.4. (42) 


Although the profiles of the bounces in FigJ^ and Eig. |2llook 
different, the bounce actions in Eqs. ( [3^-(|3^ and Eqs7p0|)- 
(42 1 are very similar]^ In particular, Eq. (39 1 provides a lower 


^By taking into account the 3(T error bands of the most important SM pa¬ 
rameters there is about one order of magnitude uncertainty on the scale where 

ySj = 0 E). 

“^This fact can be understood as follows: the bounces in Fig. can be 
rescaled by h{r) —> ah{ar) in such a way that they almost superimpose with 
those in Fig. ^ The coiTesponding change in the bounce action, induced by 
the rescaling Rsm ^sm/‘ 3. is small, because A varies very slowly ai'ound its 
minimum. 
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Figure 2: Same as in Fig. □ for the SM tree-level potential with A{\/Rsm) a 
negative constant. 


bound on the bounce action, since it is obtained by expanding 
/l(l/7?) around its minimum. 

To conclude, the good agreement between the two proce¬ 
dures for the determination of the decay probability justihes 
the approximation made in lfT2ll of taking a constant T < 0 for 
the leading order SM bounce. This was not completely obvi¬ 
ous a priori, since computing the bounce “does not commute” 
with the running of A. Of course, the potential with running A 
in Eq. ( [35] ) only captures log-enhanced corrections to the tun¬ 
nelling rate, while a complete one-loop calculation requires the 
determination of the SM action functional around the leading 
order bounce configuration IfT^ . 


2.3. Gauge independence of the tunnelling rate 


A question which is directly related to the calculation of 
quantum corrections to the SM vacuum decay rate is that of 
gauge invariance. Indeed, if one naively takes into account loop 
corrections to the tunnelling rate by computing the bounce via 


the RG improved effective potential in Eq. (35 1 the result will 
look gauge dependent. The dependence on the gauge-hxing pa¬ 
rameters in /left (e.g. in the Eermi gauge EOlU is two-fold: it 
originates from the fixed-order expression of the effective po¬ 
tential, and from its running via the anomalous dimension of 
the field h. 

Erom this point of view, gauge dependence (even if numeri¬ 
cally small) is a good thing, since is telling us that we are com¬ 
puting something in the wrong way. As shown in Ref. lfT2l . the 
divergent corrections to the bounce action are formally gauge 
independent, being directly related to the beta function of A. A 
crucial role in order to achieve the cancellation of the gauge 
dependent parts is played by the kinetic part of the effective 
action, which is neglected altogether when dealing with the ef¬ 
fective potential only. 

More generally, the gauge independence of the tunnelling 
rate directly follows from the Nielsen identity Il2ni22l 


dS 


eff 


df 


f 


d^x- 


6S 


eft 


5h(x) 


K[h{x)l 


(43) 


where f denotes the gauge-hxing parameter and is a func¬ 
tional of h whose expression depends on the gauge hxing. 

The physical implication of Eq. (j4^ is clear: the effective 
action is gauge independent when evaluated on a conhgura- 
tion which extremizes it. Hence, the bounce action is formally 
gauge independent as well. 

The proof of the gauge independence of the tunnelling rate 
can be carried out in perturbation theory by means of a loop 
expansion of the Nielsen identity, along the lines of Ref. 12^ . 
In practice, however, the cancellation of the gauge dependent 
parts in an explicit calculation might require some care. As re¬ 
cently observed in Refs. Il24l |25]| . the usual loop expansion is 
not the consistent one for the SM, where A ~ H (as in the orig¬ 
inal Coleman-Weinberg (CW) model E^ l in order for the top- 
Yukawa corrections to destabilize the tree-level electroweak 
vacuum. Consequently, such a modified loop expansion must 
be properly taken into account in order to observe the gauge 
independence of the SM tunnelling rate in perturbation theory. 
In particular, this entails the resummation of a particular class 
of daisy diagrams which (as observed in Ref. lIZTl l is connected 
with the resummation of IR-divergent Goldstone loops ll^l29l . 

Einally, we notice that for CW-like potentials, where the ab¬ 
solute minimum is radiatively generated (as in the SM), the 
standard bounce formalism requires some modifications ESI. 
On the other hand, given the fact that for the measured values 
of the SM parameters the lifetime of the electroweak vacuum 
turns out to be much larger than the age of the Universe, preci¬ 
sion calculations of the SM tunnelling rate, although important, 
are not crucial at the moment. 

2.4. Tunnelling without barriers 

The approximation of taking A a negative constant might still 
appear rather odd, since it corresponds to a tunnelling process 
from a potential of the type Alf^, with no barriers and a maxi¬ 
mum in h = 0. It is well known, however, that the absence of 
a barrier in the scalar potential is not necessarily a problem in 
held theory, due to the presence of an extra barrier originated 
by the gradient of the bounce M- In this Section we want to 
explicitly check this statement in the case of the SM. 

The tunnelling process in held theory has to be understood as 
a transition between two spatial held conhgurations at different 
euclidean times r, and Tf 

lim h(x,T,)-v, lim h(x,TA-v, (44) 

T;—»-oo T/^+oo 

where v denotes the false vacuum. The bounce action entering 
the expression of the tunnelling probability in Eq. 0 can be 
recast in a way that resembles the analogous one in quantum 
mechanics 


S[h] = 



dr K(t) yj2U[h], 


where the factor 



(45) 


(46) 
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yields the correct normalization of the path length, and 


U[h\^ JiPx + V{h{x,T)) 


(47) 


plays the role of the potential energy as in ordinary quantum 
mechanics. 

It is an instructive exercise to verify the existence of an ac¬ 
tual barrier in the case of the SM potential with d < 0. For 
simplicity (and in order to proceed analytically) we take the 
mass parameter m = 0 in the scalar potential. Starting from 
the (9(4)-invariant bounce solution Eq. ([T3]), a straightforward 
calculation yields 


T[h\ = J d^x^(Vh(x,T)) 


2 2ti^ 


Vl + ^2/ 






\^ + w 


and 


Kir)^ 


2n 




(48) 

(49) 

(50) 


Following then the definitions in Eq. ( |47] i and Eq. ( |45l l, we fi¬ 
nally get 




1 + 


1 - 


1 + F 


(51) 


and 

S[h]^ 2 drKir) yl2U[h] = —, (52) 

Jo 5|/1| 

which reproduces the correct result for the SM bounce action. 

The three quantities T[h], V[h] and U[h] are plotted in Eig.|^ 
U[h] as a function of t can be interpreted as the potential energy 



Figure 3: Shapes of the gradient {T[h]) and potential (V[/?]) barriers and their 
sum {U[h]) as a function of 

along the path which minimizes the euclidean action. Thanks 
to the positive gradient contribution T[h] we see that there is 
a barrier even for A constant and negative. Notice the correct 


asymptotic behavior of U[h], which tends to zero as t ^ oo 
since h is approaching the false vacuum v ^ 0. On the other 
hand, U[h] =0 for t - 0 corresponds to complete barrier pen¬ 
etration. Since the point t = 0 coincides with t - 0 {t - -It) 
the bounce solution can be analytically continued in Minkowski 
space, so that the system evolves towards the true minimum fol¬ 
lowing the classical equation of motion. 

Erom the above discussion one can also draw another im¬ 
portant conclusion: the largest energy scale relevant for barrier 
penetration is not that for which V = 0 (i.e. the instability scale 
of the SM effective potential), but rather the value of the bounce 
in its center h{T - 0) which corresponds to zero potential en¬ 
ergy U. 


2.5. Violation of scale invariance by mass terms 

In the previous sections (as in most of the existing literature 
on this subject) the lifetime of the metastable vacuum of the SM 
was computed neglecting the mass term of the Higgs boson, on 
the basis that the electroweak scale, of order 10^ GeV, is much 
smaller than l/R, where R is the typical size of the relevant 
bounce. We now wish to discuss this approximation in some 
detail. 

Let us consider the action 


S[h] = 



-dnhduh H— m^h^ H— Ali^ 

2 ^ ^ 2 4 


(53) 


where nf > 0 and d < 0. It was pointed out long ago Il3n that a 
bounce in this case does not exist. The easiest way to see this is 
to perform a scale transformation, defined in Eq. 0. We get 

5 [/z] ^ 5 [ha\ + ^ J d^xh^ix), (54) 

which cannot be stationary upon scale transformations unless 
h^O: 

dS[hJ 2 r A .2, , 

—-- = -m \ d xh (x). (55) 

da a=i J 

This phenomenon is well known in the context of studies of 
instanton gauge field configurations. 

Nevertheless, it is reasonable to think that, even in the pres¬ 
ence of a mass term, instanton configurations of the scalar fields 
should exist, provided they are characterized by a length scale 
R such that m « 1 /R. Roughly speaking, such a solution of the 
field equation is expected to be a function of r, approximately 
constant for r < R and approximately zero for r > R. Fur¬ 
thermore, if no other mass scale is available, the value of the 
bounce for r < R, h(Q), is proportional to 1 /R for dimensional 
reasons. Hence, the characteristic scale of h may be identified 
by the spacetime integral of a local operator, function of h. For 
example. 


f 


d\lfix) ~ h"(0) 


f 


r^dr- 


R^ 


(56) 


Based on these intuitive considerations, it is suggested in 
Ref. II 3 TI to perform a minimization of the action functional 
in which the minimum configuration is constrained to be char¬ 
acterized by a scale R much smaller than 1 jm. The constraint is 
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introduced by means of a suitable Lagrange multiplier cr, i.e. by 
adding to the action a term 



d\Oih)-cR^-” 


(57) 


where 0(h) is a local operator of mass dimension n A, for 
example 0(h) - h^, and c a constant. It is shown that with 
this modihcation an instanton appears, called a constrained in- 
stanton, which differs from the bounce of the massless theory 
by powers of m^R^, times possibly logs of inR. Thus, the cor¬ 
rections to the usually adopted approximation are indeed small, 
provided inR « 1. 

The mechanism which restores the existence of a bounce in 
the massive theory is illustrated in the Appendix of Ref 11311 
for the choice 0(h) = h^. The key point is that the constraint 
Eq. (57 1 , with cr > 0, has the effect of generating an abso¬ 
lute minimum (the true vacuum) of the scalar potential, which 
would be unbounded from below with cr = 0 and A < 0. Ex¬ 
plicitly, the new scalar potential 


Vc(h) - -m^h^ -t- -h^ + crh^ 


(58) 


3. Violation of scale invariance by new interactions 


We now turn to a discussion of the effects of SM extensions 
characterized by the appearance of an explicit new energy scale 
in the effective potential. Eollowing Refs. 0121111 , we describe 
the effect of generic NP occurring at high-energy scales by in¬ 
troducing two non-renormalizable terms in the scalar potential. 


1 

V(h) ^ -Ah 
4 


4 


—i 

6M2 



(63) 


Assuming the new couplings (/Ig g) to be of order one, the extra 
terms affect the potential only for held values of order M. The 
scalar potential in Eq. ( [6^ should be understood just as a toy 
model with no claim of being realistic. It is however sufficient 
in order to highlight some basic features of the tunnelling rate 
in the presence of new interactions. A more realistic expression 
of the modihed potential in presence of NP is obtained by con¬ 
sidering a full tower of higher-dimensional operators, namely 


- 2,2 


V(h) = -A2M^h 


-Ah^ 

4 


z 

n=3 


42n 


u2n 


2nM^" 


(64) 


with A < 0, > 0 has a local minimum at h - 0, with y(0) = 

^2" 


0, and an absolute minimum at h (for m‘^cr <si 1). The 

presence of the constraining term locally restores the scale in¬ 
variance of the action: 


da 


(S[ha]+SAha]) 


a=\ 


■f 

f 


d^xh^(x) + 2cr 


f 


A 0 2cc^ 


d‘^xh^(x) 


(59) 


which is zero for 

m^R^ 


2c 


I 


d^x h^(x) 


(mRfR^ R^. 


(60) 


This issue is not directly relevant in the SM; indeed, as we 


have seen in Sect. 2.2 scaling violation induced by radiative 


corrections have the effect of selecting a bounce of size Rsm 
10'*^ GeV"'. Explicitly, with 


V(h) = ^A(h)h"^ 


(61) 


we have 


1-5 [hj 

oa 


a=l 


-I! 


- d'^xMh)h(x), 


(62) 


and is zero around yu = 1 /Rsm ~ 10*’ GeV. As a con¬ 
sequence, scale invariance is locally restored, and a bounce of 
size R ~ Rsm is found. The effects of quadratic and cubic terms 
in the potential are suppressed by powers of vRsm ~ 10 *^ and 
can be safely neglected. Nevertheless, the simple example dis¬ 
cussed above shows that the explicit violation of scale invari¬ 
ance (for example, by the introduction of non-renormalizable 
terms in the scalar potential) may lead to complications in the 
calculation of the vacuum tunnelling rate. 


where we have also inserted a mass term due to the presence 
of a physical threshold at the scale assuming that, by a 
suitable hne-tuning, the electroweak scale is kept at the correct 
phenomenological value. 

Most of the following discussion and results apply to the gen¬ 
eral potential in Eq. ( |64l i, the key point being the explicit break¬ 
ing of scale invariance characterized by the scale M|^In princi¬ 
ple, effective operators involving derivatives of the Higgs held 
should also be included, since they inhuence the determination 
of the bounce as well. By employing a suitable held redehni- 
tion, it can be shown fi2\ that they can affect the bounce only 
when their dimension d exceeds d - 6. 

In Refs. EH El the mass scale M, required by the presence 
of non-renormalizable terms, is taken to be of the order of the 
Planck mass. Here we postpone for a while the choice of a def¬ 
inite value for M, keeping the discussion on general terms. We 
will however be especially interested in the case M » 1 /Rsm, 
since for M < 1 /Rsm the SM bounce gets trivially modihed and 
the lifetime of the electroweak vacuum can be any with respect 
to the SM case. 

Because of the dependence on M of the effective potential, 
the action is no longer scale invariant. Hence, the existence of a 
bounce for generic choices of the new couplings A 2 „ is no longer 
guaranteed. Eor example, with all A 2 n — 0 except A^ the action 
can only be stationary upon scale transformations if h = 0. We 
have indeed 

S[h]^S[h,]^S[h] + ^ Jd\h\x), (65) 


^We thank Ulrich Ellwanger for pointing out this effect to us. 

^The whole discussion applies as well to the case where such effective in¬ 
teractions arise within an explicit renormalizable model with new degrees of 
freedom, such as the one proposed in Ref. (9)- 
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which implies 


dS[h,] 


da 


a=l 


^6 


/ 


(fxh^ix). 


( 66 ) 


However, in analogy with the case of a mass term discussed 
in Section |2.5| one may argue that a constrained instanton 
with R » 1/M may still exist. Furthermore, if more non- 
renormalizable couplings are different from zero (e.g. both 
A(, and dg 0), the r.h.s. of Eq. ( [66| contains more terms. In 
such case one can conceive the possibility that these terms com¬ 
pensate each other (assuming some of the couplings have op¬ 
posite signs) for instanton solutions characterized by ~ 1 /M. 
On general grounds, we thus expect two classes of instanton so¬ 
lutions (when they exist); those characterized by » 1 /M and 
those characterized by ~ 1 /M. 

With a generic potential of the type of Eq. ( |64| ), the field 
equation cannot be integrated analytically, and one has to adopt 
numerical methods (which we are going to do below). How¬ 
ever, a somewhat deeper understanding of the tunnelling pro¬ 
cess in the presence of NP can be achieved by making contact 
with the case of the pure SM. To this purpose, it is useful to 
rewrite the modihed potential Eq. ( |64| i in terms of an effective 
quartic coupling]^ 


v(h) = 


with 


Aeff(h) - 


2A2M^ 

h?- 


■ A + 


n=3 


A2n 


2nM^” ^ 


.2n-4 


(67) 


( 68 ) 


where, unless otherwise specihed, we take the quartic coupling 
A to be at its minimal negative value (as in the SM case) and 
we neglect the logarithmic running for the NP couplings, since 
their are anyway unknown. Note that the coupling A 2 becomes 
active only above threshold h> M. 

If we were able to argue that the dominant contribution to 
the bounce action is provided by taking the argument of defr 
at some hxed value, then the action would be approximately 


given by the SM expression Eq. (26 1 , with A replaced by the 
appropriate value of A^ff. 

It can be seen that this is indeed the case by adopting the 


approximation mentioned in Section 2.5 namely, to consider 


the bounce as a finite-action, (9(4)-invariant solution of the eu¬ 
clidean field equation /i(r), approximately constant for r < R 
and approximately zero for r > R. Let us first test this approxi¬ 
mation in the case of the SM with no mass term, where the ex¬ 
act result is known. In this case, the constant value h{Q) of the 
bounce for 0 < r < R must be proportional to l/(R VW)- The 
factor l/R arises for dimensional reasons. The factor 1/ VW 
can be understood by writing the action as 


S sm[/*] - J" 


d^x 


-d^hduh — Ah^ 
2 ^ 4 




(69) 


tial. 


'^Note that this parametrization is not directly related to the effective poten- 


where the bounce 


H{x) = ^4Ah{x) 


(70) 


does not depend on A. 

Within this approximation we get 


S SM[h] 


^27^ 

Jo 


drr 


V(h) - -h(r)V'(h) 


jjR^ 


2 

n^R* 


8 


v(hm - -h(o)v'(hm 


-\A\h\0) 


(71) 


where a is a constant, and we have used the field equation in 
the first equality. The correct value SsmW = 87r2/(3|/l|) is re¬ 
produced for = 8/ V3 =» (2.15)^. 

We now turn to the case of the modified potential Eq. 

We observe that, within the same approximation that led us to 
Eq. (1^, the effect of a scale transformation Eq. ( [TT] i on the 
modified action is 


dS[h,] 


da 


= r 

a=l J 


dA^ffih) , 
4 


^ \ ~ h{0)A'^g(h(0)). (72) 


dh 


Hence, in order for the action to be stationary under scale trans¬ 
formations, !i{Q) should be chosen so that 


4' (h(0)) = 0 


(73) 


locally restores scale invariance in a neighbourhood of h(Q) (the 
case h{Q) = 0 is obviously uninteresting). If such a value of h(Q) 
exists, then 


S [h] ^2 k 
tJR^ 


r 


drr 


V(h) - ^-h(r)V'(h) 


2 

n^R^ 


V(h(0)) - -h(0)V'(hm 


1 


A,s(hmh\o) - 


8 


2d4 


tt'^R 


-A,s(h(Q))h\Q), 


(74) 


which is precisely what we would get in the SM with no mass 
term, and A replaced by /leff(/i(0)), see the third line in Eq. ( |7T] i. 
We conclude that the leading contribution to the bounce action 
is given by 

_ _ (75) 


3|4efF(/!(0))|’ 


provided 


h(0) = 


R V|4eff(/l(0))| 



— . (76) 


As an example, let us consider the potential in Eq. (63 1 . We 
have in this case 


2 I // 

/IpfrC/z) — A + ~Af, -r- + —/Is- 7 , 

^ '' 3 ^2 2 M4’ 

, 2 h I h^\ 


(77) 

(78) 
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We can distinguish three basic scenarios. 


I. One possibility is that both dg and dg are non-negative and 
at least one of them is different from zero. In such case the 
potential is bounded from below and, for any value of h, 

detf(/!) > A{h). 


In the limit where we assume d = -|/l| constant, is 

always different from zero; local scale invariance is hope¬ 
lessly lost, and a bounce cannot be found. These expecta¬ 
tions are conhrmed by our numerical analysis, described 
in Appendix A Actually one hnds that a bounce with a 
SM-like action exists for hnite values of JCmax, but its size 
R grows to inhnity as Xmax is sent to inhnity. Since h{Q) is 
inversely proportional to R, it follows that the bounce is in 
fact zero everywhere. 


However, the absence of a bounce turns out to be an arte¬ 
fact of choosing d constant in Eq. ([68]l. With d running 


there is an extra compensating contribution in Eq. (78 i, 
and indeed a bounce of hnite size with a SM-like bounce 
action can be found numerically, provided M » 1 //?sm ~ 
lO'^ GeV, the scale at which d has a minimum. Hence we 
conclude that in such a case the decay probability of the 
false vacuum is not modihed with respect to the SM value. 


II. A second possibility is that dg and dg have opposite signs, 
but the potential remains bounded from below, hence dg < 
0 and dg > 0. In this case, the two non-renormalizable 
terms compensate each other and restore scale invariance 
locally (in the neighborhood of held conhgurations char¬ 
acterized the scale 1/M). Indeed, is zero for 


h = h{Q) = M 



(79) 


and 

4ftW0)) = -|d|-^. (80) 

ydg 

The instanton size is now immediately read off Eq. ( [7^ : 

a a \ 1 

_ _ ___ 

“ hiO) V|deff(h(0))| “ M ~ M' 

\ S/lg y ' 9/^8 / 

( 81 ) 

Two important points are to be noted. Eirst, we have in 
this case 

|deff(/r(0))| = |d| + ^ > |d|, (82) 

ydg 

and therefore 

S[h]<SsM[h]. (83) 

As a consequence, the decay probability of the false vac¬ 
uum can only be increased with respect to the pure SM. 
Second, the dominant contribution to the tunnelling rate is 
independent of the scale M of NP, because the action is. 
This conclusion has been reached in the context of an ap¬ 
proximate calculation, but it is easy to see that it holds true 


A 

dg 

dg 

h 

R 

S{h\ 

const 

<0 

< 0 

0 

CO 

0 

const 

<0 

> 0 

+ 0 

~ 1/M 

< 5sm[/i] 

const 

>0 

< 0 

0 

CO 

0 

const 

>0 

> 0 

0 

CO 

0 

run 

<0 

< 0 

0 

CO 

0 

run 

<0 

> 0 

+ 0 

~ 1/M 

< 5sm[/i] 

run 

>0 

< 0 

0 

CO 

0 

run 

>0 

> 0 

+ 0 

Rsu 

S sm[/i] 


Table 1: Summary of the bounce solutions associated with the potential in 
Eq. ^631 Only the case M a> 1 /Rsm is considered, while for M < 1 /Rsm the 
bounce action can be any with respect to the SM one. 


in general. Indeed, one may dehne dimensionless coordi¬ 
nates and held by 

h(r) 

X = Mr, h(x) = (84) 

so that 


V{h) = 


1 


-Ah^ + 


Z 

«=3 


J 2n 

2n 


= M^Vih). 


(85) 


Then 


S[h\^2y drr^ 
Jo 


r 


V(h) - -h(r)V'(h) 


= 27r‘ 


r 


dxx~ 


v(h) - -kxw'iJi) 


( 86 ) 


which is independent of M. This is an important point; 


whatever modihcation may the potential in Eq. (64 1 induce 


on the tunnelling rate, this is independent of the scale M 
up to subleading pre-exponential corrections. 


III. The last case to be considered is the one where the poten¬ 
tial is unbounded from below at large held values. When 
both Tg and Tg are negative (or at least one is negative and 
the other is zero), it is clear that the bounce does not exist 
because of the lost of scale invariance. 

The case dg > 0 and dg < 0 is more subtle; here a non-zero 
solution of d'jj:(/r) = 0 does exist, but it corresponds to a 
maximum of deft and not to a bounce conhguration. Actu¬ 
ally the fast drop of the potential at large held values makes 
the whole problem of hnding a bounce solution ill-dehned 
in this case; there exist “rolling solutions” that destabilize 
the electroweak vacuum characterized by a much shorter 
time scale M- 

These results are summarized in Table [T] with either constant 
or running dj^ To complete this study, we show in Eig. the 
comparison between the numerical determination of the bounce 
action and the analytical approximation in Eq. ( [75| ). The agree¬ 
ment is generally quite good. Hence, one can take advantage 


similar study of the tunnelling rate as a function of dg and dg was per¬ 
formed in Ref. (H 
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A=-0.01473, A8=2.1 



Figure 4: Numerical vs. analytical (approximated) determination of the bounce 
action as a function of A(,. The relative difference is always within 30%. For¬ 
mally, for A constant and > 0 a bounce does not exist; the result in the plot 
for /l6 > 0 has to be understood for A running. 



Figure 5: Evolution of /lefr(/2(0)) (including the running of d) for three different 
values of the NP scale, M = 10^ GeV, Mp and and for three different 

choices of the NP parameters /l6,8- The instability bound is defined by the 
inequality |/lcff(/i(0))| > ^ logfrUco)) ’ = 4-35 X lO'^ sec is the age 

of the Universe. 

of the /leff language in order to describe in an essential way the 
tunnelling process, as exemplified in Fig. where some po¬ 
tentials corresponding to cases I and II are considered. From 
that it is clear that modifications of the SM potential above 
1/^SM ~ 10*^ GeV may only shorten the lifetime of the elec- 
troweak vacuum. 

The above discussion can be generalized to the case of more 
non-renormalizable terms, including derivative operators, the 
three main categories being defined by the stability (I and II) or 
instability (III) of the potential at large field values, and by the 
presence of a new minimum around the scale M (case II). 

Summarizing, the explicit breaking of scale invariance at 
energy scales above 1//?sm ~ 10*^ GeV may lead to a new 
bounce, whose characteristic scale is necessarily related to that 
of NP (case II). However, this happens only if the NP modifies 
the ground state of the theory, generating a new deep minimum 


around the new scale M. In such case, the effect of NP is that of 
opening a new decay channel for the electroweak vacuum and 
the tunnelling rate can only increase with respect to the SM cal¬ 
culation. It is then obvious that in this case one is not analyzing 
anymore the stability of the SM, but that of a different theory, 
with a completely different ground state. 

The above argument also explains why there is an apparent 
violation of the decoupling theorem, according to which one 
expects new-degrees of freedom characterized by a scale M » 
1 /Rsm to be irrelevant for the analysis of the vacuum stability. 
This apparent paradox is due to the fact that the modifications 
of the effective potential introduced by means of Eq. ( [6^ are 
not only related to the appearance of high-frequency dynamical 
modes (for which the decoupling still applies), but they also 
imply a drastic modification of the ground state of the theory 
(invalidating the decoupling argument). 

We finally note that the modified effective potential of the 
type in Eq. ( [6^ , leading to a fast tunnelling rate (case II), is not 
a well motivated UV completion of the SM potential close to 
the Planck scale. On the one hand, except for fine-tuned values 
of Te and dg, the lifetime of the electroweak vacuum turns out 
to be extremely fast, in sharp contradiction with the existence of 
the present Universe. This implies that such modified potential 
cannot be considered as a phenomenologically viable UV com¬ 
pletion of the SM potential. On the other hand, if M is close 
to the Planck mass, truly gravitational effects 041 cannot be ig¬ 
nored. As recently pointed out in Ref. Osl . the latter suppress 
the tunnelling rate by many orders of magnitude]^ 

4. Conclusions 

The decay probability of the electroweak vacuum under 
quantum tunnelling is intimately related to the breaking of scale 
invariance in the Higgs effective potential. The latter can be 
characterized by an effective h* interaction, both within and be¬ 
yond the SM, which is not scale invariant beyond the tree level. 
To a good approximation, the size of the bounce and the tun¬ 
nelling probability are determined by the energy scale at which 
this effective coupling, Aggih), reaches its minimum value and 
by its value at the minimum, respectively. In the absence of 
NP, the breaking of scale invariance occurs dominantly via SM 
radiative corrections. The latter selects a leading bounce char¬ 
acterized by the scale 1/Rsm ~ 10'^ GeV, that implies a suf¬ 
ficiently long lifetime of the electroweak vacuum compared to 
the age of the Universe. 

New degrees of freedom at high energies can, in principle, 
introduce a new explicit breaking of scale invariance in the ef¬ 
fective potential. This, in turn, can shift the energy scale where 
Aeffih) reaches its minimum value and, most importantly, can 
modify such minimum value. If the new degrees of freedom ap¬ 
pear above the scale 1 /Rsm, we can distinguish two cases: those 
where NP simply stabilizes the SM potential, and those where 


^Gravitational effects in minimal Einstein gravity tend to slow down the 
tunnelling rate (HI However, as long as \/Rsm ^ these corrections are 
small and hence the presence of gravity does not inficiate the results obtained 
within the SM in isolation (32) 
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NP introduces new decay channels for the electroweak vacuum. 
In the first case, characterized by the absence of new minima 
for the lifetime of the electroweak vacuum remains un¬ 

changed. In the latter case it can only decrease with respect to 
the pure SM. In particular, a shorter lifetime occurs only if the 
new degrees of freedom drastically change the ground-state of 
the theory, introducing a new minimum for A^gih), such that it 
does not make sense anymore to speak about the stability prob¬ 
lem of the SM potential. 

Given these arguments, we can conclude that the sensitiv¬ 
ity of the tunnelling rate to the possible UV completions of the 
model does not invalidate the vacuum stability analyses per¬ 
formed using the pure SM potential. The scope of the latter is 
answering the following well-defined physical question: Does 
the extrapolation of the SM up to the Planck scale necessar¬ 
ily implies the existence of NP below such scale? According 
to the present experimental values of ot/, and m,, we can state 
that the answer to this question is no IT] |3- Were the top mass 
e.g. 180 GeV, the answer would have been different since NP 
would have been necessarily implied well below the scale where 
A reaches its minimum SM value, regardless of the physics in 
the deep UV which cannot improve on stability. We finally 
stress that the negative answer to the above question does not 
necessarily imply the absence of NP up to the Planck scale: it 
only says that the model is compatible with the absence of NP 
up to the Planck scale. Similarly, from this analysis we cannot 
infer that any UV completion of the theory at the Planck scale 
is compatible with the stability of the electroweak vacuum, but 
only that it is possible to build UV completions that do not con¬ 
tain new degrees of freedom below the Planck scale. 


Acknowledgments 


with boundary conditions 


lim h(x) - vjM, 

x—^oo 


dh(x) 


dx 


= 0. 


t=0 


(A.3) 

(A.4) 


A numerical solution can be found by means of the shooting 
method (see also ||8l|33l|36l), which consists in tuning the ini¬ 
tial condition aq in such a way that Eq. ( |A.3| l is satisfied at the 
boundary. In practice, since we cannot start at the singular point 
.x = 0, we Taylor-expand the solution around the origin and de¬ 
fine the Cauchy problem 

/i(xmin) = ao + ^V'(M0))xA„, (A.5) 

b (.Vjjiin) — 4^ (^(0)).^min5 (A.6) 

where we kept up to quadratic terms in Xsain- By evolving the 
solution h{x; ao) up to .x = Xaiax, we require 




(A.7) 


which determines the value of do which satisfies the boundary 
condition at x — x,^^. 

Analogously, by Taylor-expanding the solution at the first 
non-trivial order in 1 /jCmax for large Xmax, we get (m -0 case) 

^(^max) = 4^, (A.8) 

■^max 

which is useful in order to have an analytical control over the 
asymptotic solution. On the other hand, if the mass term is kept 
in the potential (m 0 case) the asymptotic solution is (see 
e.g. Appendix of llJTll ') 
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^(■^max) “ ^oo Ai(c.\:jnax) ~ ^oo 

•^max •^max 



(A.9) 


where e - m/M and in the last step we used the asymptotic 
expression of the modified Bessel function Kfx) for large x. 

The bounce action is determined by integrating the profile of 
the bounce h(x', do) between Xmin and Xmax 


27T" 


f 

*yXmi 


S[h]^ I xMx 


V(hM) - ^iMV'ihM) 


(A. 10) 


Appendix A. Numerical determination of the bounce 

For the numerical analysis it turns out to be convenient to 
work with adimensional variables 

X = Mr; h{x) = —; V(h) = (A.l) 

M M* 

where M is an arbitrary mass scale. The bounce equation reads 
Tf'ix) + -h'(x) = V'(h), (A.l) 

X 


where we applied the equation of motion and performed the 


change of variables in Eq. (A.l I. Finally, the algorithm is it¬ 
erated by choosing increasingly smaller (larger) values of Xmin 

(■^max)- 
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